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Abstract 

Range and load play keys on the problem of attacking on links in random scale-free (RSF) 
networks. In this Brief Report we obtain the relation between range and load in RSF networks 
analytically by the generating function theory, and then give an estimation about the impact of 
attacks on the ef ficiency of the network. The analytical results show that short range attacks 
are more destructive for RSF networks, and are confirmed numerically. Further our results are 
consistent with the former literature (Physical Review E 66, 065103(R) (2002)). 
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Attacking on complex networks, especially in the context of the Internet and biological 
networks, has been an interesting issue, and different aspects of attacking have been analyzed 



recently Many works focus on attacks ^n nodes, and the strategies provided include 

random attacks, degree-based attacks, etc. [2]. Also, some works consider attacks on links, 

II 

and the strategies include the range-based attacks, load-based attacks, etc. p]. 

Adison E. Motter, Takashi Nishikawa and Ying-Cheng Lai jj] studied attacks on links 
in scale-free networks basing on range. Range is introduced by Watts [s^ to characterize 
different types of links in networks: the range of a link lij connecting nodes i and j is 
defined as the length of the shortest path between the nodes i and j in the absence of kj. 
The small-world model introduced by Watts and Strogatz [6] (WS model) is more sensitive 
to attacks on long range links connecting nodes that would otherwise be separated by a 
long distance. It is not true for many scale-free networks, though most of them also have 
a short average path length like the WS model. Adison E. Motter, Takashi Nishikawa and 
Ying-Cheng Lai found that short-range links other than long-range ones are vital for efficient 
communication between nodes in these networks. They argued that the average shortest 
path is a global quantity which is mainly determined by links with large load, where the 



3, 



And 



load of a link is defined as the number of shortest paths passing through this link 
for scale-free networks, with exponent in a finite interval around 3, due to the heterogeneous 
degree distribution, the load is on average larger for links with shorter range, making the 
short-range attacks more destructive. 

In this brief report, employing the generating function theory, we first derive an approxi- 
mate relation between R{ki, ^2) and L{ki, ^2) for RSF networks analytically, where R{ki, 
and L{ki, ^2) are defined as the expected value of range and load respectively for links be- 
tween nodes with given degree ki and ^2. We then give an estimation about the decrement 
of ef ficiency as a function of R{ki, ^2) and L{ki, /C2), showing that short-range attacks are 
more destructive for RSF networks. Numerical simulations are also performed, confirming 
our analytical results. 

To study ranged-based attacks on links in RSF networks, we measure the ef ficiency of 
he network as each link is removed. The ef ficiency of a network with size N is defined as 

a: 



^ N{N 

where dij denotes the length of the shortest path between the node-pair {i,i), the sum is 



over all pairs of nodes in the network. The ef ficiency defined above has a finite value even 
for disconnected networks, and larger values of E correspond to more efficient networks. 

When a link is removed from the network, the ef ficiency of the network generally de- 
creases. The decrement of ef ficiency involves two quantities: (1) the number of node-pairs 
whose geodesic lengths increase; (2) the average increment of the geodesic lengths of these 
node-pairs. The first quantity is related to the load of the removed link, and the second 
quantity is related immediately to the range of the removed link. 

For RSF networks, the expected value of the geodesic length of node-pairs with given 
degree ki and k2 is 

d{kik2) = ^^P* (^1,^2) , 

i=l 

where p^{ki k2) is the probability that the node-pair with given degree ki and k2 has a 
geodesic length i. For RSF networks, we have [lO] 

At,M - ^. (2) 

where is the number of nodes in the network, zi is the average number of first neighbors. 
By the generating function formalism, we can obtain [1 1|] , 

\n{z2/zi} 

where Z2 is the the average number of second neighbors. Accordingly the expected diameter 
of RSF networks is HI] 



ln(iVAi) 

Since the RSF network is totally random in all aspects other than the degree distribution, 
R{ki,k2) is thus equal to the expected value of geodesic length of nonadjacent node-pairs 
with given degree ki — 1 and /c2 — 1, that is 

P^iki-l,k2-l] 

■I — 
1 

i=2 



I.e., 



. . djki - 1, fc2 - 1) - p\ki - 1, fc2 - 1) . . 

1 -pl(/Ci - l,fc2 - 1) 



Combining Eq. ([2]), ([3]) and (j5]), we can obtain 



ln(jVzi/((fci-l)(fc2-l))) 
2Nzi 



Furthermore, we assume that the network is spare, and can be seen as a tree with expected 
diameter D. Consider a hnk Uj connecting node i and j, where i has a degree ki and j has 
a degree /c2- When removing /j^, the network can be regarded as a tree Tj rooted as i or Tj 
rooted as j, both of which have a depth oi D — 1. Staring from the root z, the first layer 
has ki — 1 nodes, the second layer has Zi{ki — 1), and the mth (0 < m < D) layer has 
z^~^{ki — 1) nodes. Similarly, the mth layer of Tj has z'^~^(k2 — 1) nodes. The geodesic 
path from nodes in the dith [di < D — 1) layer in Tj to nodes in the d2th {d2 <= D — 1 — di) 
layer in Tj in are expected to pass kj, which has contribution of 1 to the load of kj. Thus 
the expected value of the load of kj is 



L{k,,k2) = - 1) ^'' ^^1 ^ + 1)(^2 - i){zi - 1)""-'-' 



+m - 1) ^'^"^^''''"^ + 1) + ih - - 1)^-^ 

= ( ^2 iz,-2y )ik,-l)ik2-l)) 

+{k, + k2- ^ + 1. (7) 

Zi — z 

When ki » Zi,k2 » zi, the above equation can be rewritten as 

L{ku h) = _ \ki - l)(/^2 - 1), (8) 

showing that load is directly proportional to the product of {k,\ — 1) and {k2 — 1) when fci 
and ^2 are large enough. For simplicity, we rewrite Eq. (JHl) as 

L{k^M)=c{kr-l){k2-\), (9) 

where c is the coefficient '^^^^'^^^^^^f>^^]!^^^'^^^^ ±1. 

The above analytical results can be numerically verified in the following. In our example, 
we take = 10000, A = 3.5, k^ = 6. We plot R{ki, ^2) in Fig.l, and L{ki, ^2) in the inset 
of Fig.l. From the inset of Fig.l, it can be seen that load is directly proportional to the 
product {ki — l){k2 — 1) when {ki — l){k2 — 1) are large enough. 

Combining Eq. ([6]) and ([9]), we can see the load and the range have a negative correlation, 
that is 

ln{cN-zi/L) 

R(^L) ^ 1 + J^i£lZ|^. (10) 



This expression gives an estimation of the relation between R and L, suggesting that short 
range hnks are expected to be passed through a large number of shortest paths. Numerical 
verification is presented in Fig. 2, where theoretical estimation and numerical results are 
plotted. When a link lij is removed from the network, the decrement of the ef ficiency of 
the network is approximately: 



AE 



2 >^ R{t,j)-l _ 2iRit,j)-l)L{t,j) 



where the set F is all the node-pairs whose shortest length should increase as a result of the 
removal of link lij, R{i,j) and L{i,j) are the range and load respectively of lij. Thus the 
product {R — 1)L is a natural quantity to characterize the impact of removing a link on the 
ef ficiency. For RSF networks, 

lnicN-zi/L{i,j)) 

AE^hL{^,J)^^^^^, (12) 

^ 2Nczi 

where h is ^2jv(Af-i) • is an increasing function of L, thus a decreasing function of R. It 
can be concluded that links with small range are more important for the efficiency of RSF 
networks. 

In summary, by investigating the expected range and load of links in RSF networks, 
we obtain an approximate analytical relation between range and load, and then give an 
estimation of the impact of removal of links on the ef ficiency . Thus we prove analytically 
that attacks on short-range links are more destructive for RSF networks. An insufficiency in 
our work is that R{L) has a significant error comparing with numerical results. R{L) given 
in Eq. (fTOl) is highly sensitive to some quantities such as c and Z2. These quantities are not 
obtained accurately due to the approximations used in our analysis, such as the tree-like 
approximation necessary for the generating function method. Anyhow, the analytical results 
in this brief report give a reasonable description for the trend of the true relation between 
R and L for RSF networks. 
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Fig.l Average range as a function of the product {ki — l){k2 — 1) in RSF networks 
with N — 10^, A = 3.5, mo = 6, and the maximal degree mmax = 500. The sohd hne is 
the theoretical curve and the hollow squares are simulation results. Inset: Average load 
(represented by B) as a function of the product {ki — l)(/c2 — !)• Numerical data are 
obtained from 100 realizations. 

Fig. 2 Average load as a function of range. Square: theoretical value; circle: averaged 
simulation value over 100 reahzations; (All the parameters are the same as in Fig.l) cross: 
simulation value for A = 3.0, — 5000 from Adison E. motter, Takashi Nishikawa and 
Ying-Cheng Lai's work (Physical Review E 66, 065103(R) (2002)). 
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